Measurement of Coulomb drag between Anderson insulators 
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We report observations of the Coulomb drag effect between two effectively 2-d insulating a- 
Sii-ajNba; films. We find that there only exist a limited range of experimental parameters over 
which we can measure a sizable linear-response transresistivity (p<j). The temperature dependence 
of p c \ is consistent with the layers being Efros-Shklovskii Anderson insulators provided that a 3-d 
density of states and a localization length smaller than that obtained from the DC layer-conductivity 
are assumed. 



Materials such as a-Sii-^Nba; which exhibit a disorder 
driven T — (Quantum Critical) Metal-Insulator Tran- 
sition (MIT) [H, i] have presented many challenges to 
condensed matter physics: in particular understanding 
the role of long ranged electron-electron interactions in 
the insulating phases Q and in the vicinity of the MIT 
0. Since the interplay between disorder and electron- 
electron interactions in such systems will determine the 
dominant transport mechanism, the development of ex- 
perimental techniques to separately measure these is use- 
ful. In this letter we show that the Coulomb drag effect 
allows us to directly study long ranged electron-electron 
interactions in insulating a-Sii-^Nb-r thin films. We find 
that although linear-response Coulomb drag is only ob- 
servable over a limited range of sample parameters, when 
obtainable, it offers unequivocal distinction to be made 
between alternative models for the electronic transport 
in such systems. 

Coulomb drag Q arises from the Coulomb scattering 
of charge carriers in spatially separated layers, in the ab- 
sence of charge transfer between the layers. Experimen- 
tally, the Coulomb drag effect between two layers (layer- 
I and 2) can be observed by measuring the electric-field 
(E 2 ) created in one, open circuited, layer due to a parallel 
applied current-density (ji) in the other. The (longitudi- 
nal) transresistivity, or linear-response Coulomb drag co- 
efficient, is defined as pd = —Ez/ji] whilst the measured 
total-transresistance is the ratio of the induced voltage 
in layer-2 to the applied current in layer-f, i.e. —V2/I1. 
Theoretical analyses agree that the linear-response tran- 
sresistivity between two identical 2-d layers is given by, 
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where (3 = (fceT) -1 , T is the temperature, n is the carrier 
density in a layer, x( w ? q) is the density-density response 
function of a layer, and U(u>,q) is the screened interaction 
potential between the layers. Imx(w, q) may be obtained 
from the finite wavevector conductivity via: 

lm X (w,q)=q 2 (e 2 u ) )- 1 a{u>,q). (2) 

The motivation for this study comes from predictions @ 
that the temperature dependence of the transresistance 



can be used to differentiate insulating states. In partic- 
ular, for the case of (2-d) Mott Anderson insulator bi- 
layers the low-T transresistivity should vary as pd oc T 2 , 
whilst for (2-d) Efros-Shklovskii(ES) Anderson insulators 
it should diverge when T -> as p d oc T 3 exp[(r /T) 1 / 2 ], 
where To is the ES characteristic temperature given by 
ksTo « 2e 2 /k£ and n is the dielectric constant of the 
layers. This opposite behavior of the transresistivity as 
T — > (diverging for ES- Anderson insulator bilayers, and 
vanishing for Mott Anderson insulator bilayers) allows for 
a more transparent distinction between the two insulat- 
ing states over a narrow T range than intra-layer trans- 
port measurements do — where one would obtain different 
functional T or u> dependencies with the same trend. 

In our study, samples consisted of two U-shaped 
200A thick insulating a-Sii-^Nb^ layers, separated by 
a SiO barrier (see upper left inset in figure [2]). All layers 
were fabricated using standard RF magnetron sputtering 
techniques, with the a-Sii-^Nbz layers being deposited 
using the co-sputtering technique with a rotating sample- 
holder outlined in Ref[3(. Samples were all grown on 
polished glass slides in an inert (argon) environment at 
ambient temperature and a pressure of f=a 0.50kPa. Prior 
to the layers being deposited, 600A-thick silver connec- 
tions and a silver bridge, that would connect the two lay- 
ers, were sputtered. The latter would prevent the barrier 
from breaking due to electrostatic discharge/breakdown 
during handling, characterization and cooling. Connec- 
tions to the signal generating and measuring apparatus 
were made with sputtered silver layers on the arms of the 
U-shaped layers several millimeters from the interacting- 
rcgion. To prevent oxidation the top layer was covered 
with a ?s50A-thick SiO film prior to atmospheric expo- 
sure. Thicknesses were inferred from the sputtering rate 
(measured using a quartz crystal thickness monitor) , and 
were within w 5% of those obtained from optical inter- 
ferometry measurements on samples of similar composi- 
tions and thicknesses. Estimates of the uncertainty in 
the thickness due to fluctuations in the deposition rate 
and uncertainties in the exposure time are < ±1 OA and 
< ±15A for the case of the a-Sii-^Nb^ and the barrier 
layers respectively. The barrier layer was deposited in 
discrete stages to reduce the formation of pinholes. High 
sputtering powers and brief atmospheric exposure be- 
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FIG. 1: Plot of the low-T transresistance per square 
(fl^drag /square) for various samples, at a driving-current of 
Ii — InA. (The right-vertical axis shows the ratio of the total 
measured voltage to the total driving current). Samples 1A, 
IB and 1C have an average Nb concentration of x = 0.070 
and layer separations of 200A, 100A and 55A respectively. 
Samples 2C and 3C have the same geometry as sample 1C 
but with x » 0.076 and » 0.080. 



tween stages were found to increase the barrier strength 
and durability. We were thereby able to fabricate barrier 
layers with resistances several orders of magnitude larger 
than those of the layers. 

Samples were cooled to T « 1.2K using a standard 
liquid helium cryostat. The low-T (<20K) resistance 
of the barrier layer was obtained by measuring its con- 
ductance, i.e. by applying a small DC potential differ- 
ence between terminals on either layer, and then measur- 
ing the induced current. Measurements were performed 
from several terminals on each layer and with various 
terminals and connections grounded to rule out any ef- 
fects from ground- leakage and -loops. Single layer re- 
sistance measurements were performed using standard 
4- wire techniques, for which it was assured that driv- 
ing currents were small enough to give linear-response 
coefficients. The indifference of the layer resistances to 
the grounding of terminals in the second layer provided 
further verification that inter-layer leakage and tunnel- 
ing was not significant. Transresistance measurements 
were performed using a quasi-AC technique in which 
a DC current source was programmed to flip polarity 
every ~1 second. Since in linear-response theory pd 
remains the same on interchanging the layers (even if 
the layers have different resistivities), the condition that 
(-Rdrag/square) 2 —i = (-Rdra g /square)i^ 2 as the driving- 
current is decreased, was used as a test for the linear- 
response regime. The transresistance at driving currents 
of J = InA for samples with average Nb concentra- 
tions per layer [3] of x = 0.07 — > 0.08 and layer sepa- 
rations of 50 -> 200A were found to, on average, increase 
with decreasing temperature between T ~ 20 — > 3K. 
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FIG. 2: Resistance per square of layers 1 (triangles) and 2 
(circles) in sample IB. Open symbols are plotted against the 



T~ ' axis and closed symbols against the T~ ' axis. Upper 
inset: sketch of sample geometry. Lower inset: resistances 
of layers-1 and -2 (triangles and circles) compared to lower 
bounds of the barrier resistance (squares). The dashed hori- 
zontal line represents the maximum resistance that could be 
measured using Vmduced(^appiied) techniques with existing ap- 
paratus. 



At lower temperatures the transresistance would satu- 
rate or decrease in magnitude. This is shown in fig- 
ure [D for a selection of samples. In what follows we will 
present data for sample IB, where the transresistance 
entered the linear-response regime for I < InA. Due 
to the multitude of competing size effects (see below), 
it was unclear, based on the temperature dependence 
of the DC resistance alone, whether the layers are bet- 
ter described by the ES Variable Range Hopping(VRH) 
model Q, which predicts: \n(a<ic) °c T -1 / 2 , or the Mott 
VRH model 9], which, for effectively 2-d films, predicts: 
m (°dc) oc T^ 1 / 3 , at low temperatures. This can be seen 
in figure [H where we present T~ 1//2 and T -1 / 3 plots of 
the layer-resistances of sample IB. In the temperature 
range 4K < T < 15K, ES VRH with a characteris- 
tic temperature of T = 145(± < 9)K, 2-d Mott VRH 
with a characteristic temperature of To = 1890(±60)K, 
and 3-d Mott VRH with a characteristic temperature 
of T = 36, 000(±1000)K, all describe the T-dependence 
reasonably well. The localization length determined from 
the ES characteristic temperature is £ ~ 105 (±6) A. For 
each case the DC VRH transport in the layers should 
be effectively 2-d for T<20K — on account of the reso- 
nant hopping distance (r c ) being larger than the width 
of each layer (W). At these temperatures the resistance 
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of the barrier — obtained by measuring the tunneling cur- 
rent (see above) — is found to be approximately two or- 
ders of magnitude larger than the resistance of layers- 1 
and -2 (see lower inset in figure [2]) . 

The observed decrease of the layer-resistance and tran- 
sresistance observed at T<2.5K is (based on the tests 
outlined above and in Ref[ll|) neither due to grounding 
loops/leaks or tunneling/leakage through the barrier. It 
is unlikely that they are due to thermoelectric effects (as 
in e.g. [4[), given that all connections were far away from 
the interacting region. Furthermore, it is unlikely to be 
due to the layers undergoing a superconductor-insulator 
transition — which is known to occur in thin a-Sii-^Nb^ 
films with large Nb concentrations (x ~ 0.15 — 0.18 
fiol ]) — based on the high resistance per square of the lay- 
ers and the high temperature at which the transition ap- 
pears to occur. We will present a comprehensive study 
of this phenomenon elsewhere, and focus on the temper- 
ature range 2.5K < T < 15K in this letter. 

As the driving current is increased the transresistance 
decreases (see inset in figure [3]), and for I ks lOOnA, the 
transresistance is found to be more than 3 orders of mag- 
nitude smaller than at I = InA. This is likely due to 
larger driving currents both increasing the effective sam- 
ple temperature and producing significant non-linear re- 
sponses. For / < InA the drag voltage becomes noisier, 
but its average value changes by less than w 5% between 
1, 0.75 and 0.5nA (see figure [3]). We thus believe that 
at this point we have reached the linear-response regime. 
It is unlikely that phonon 12 1 and plasmon [l3| contri- 
butions are significant at these temperatures on account 
of the small layer separations and the strong insulating 
nature of the layers. 

In order to limit the number of fitting parameters 
in analysing these data, we find it convenient to plot 
the ratio (i?drag/square)/(i?i aycr /square) as a function 
of temperature. Doing so we find that our data is 
best described (see figure 0} by the 2-parameter equa- 
tion: (i?drag/square)= aT b (R\ aycl /square), with a = 
1.1(±0.1) x 10" 4 K- 6 and b = 2.0(±0.1). The latter pa- 
rameter deviates from the predictions of RefQ, which 
suggest that b = 3 for a bilayer system comprised of 2-d 
ES Anderson insulators. 

The observed discrepancy can be explained if the 
screening in the layers is not 2-d in the studied regime — 
i.e. if the response is dominated by a finite wave-vector 
q > W^ 1 (where W ~ 200A) for which the layers 
will be effectively 3-d. Since > W^ 1 > r^ 1 = 
£ _1 (^b/2') -1 ' 2 , at the temperatures of interest, this oc- 
curs at momentum transfers (q > r" 1 ) that dominate the 
transresistance (see Ref@). 

Repeating the calculations for the finite-u;, <? conduc- 
tivity in Ref[3] for 3-d systems, we find it takes the 
asymptotic forms: 
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FIG. 3: Temperature dependence of linear-response transre- 
sistance observed at driving currents I < InA for sample 
IB. Inset: At driving currents / > InA. The anomalous low- 
T increase of the / < InA transresistance is consistent with 
the layers being Efros-Shklovskii Anderson insulators and not 
Mott Anderson insulators 
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FIG. 4: Plot of the T-dependence of the ratio (7?dra g /square) 
/ (-Riaycr-2 /square) at I — InA for sample IB (circles), the pre- 
dicted slope for the case of two 2-d Mott- Anderson insulator 
layers (dotted line), the predicted slope for the case of two 
2-d Efros-Shklovskii Anderson insulator layers (dash-dotted 
line), and the predicted slope for two effectively 3-d Efros- 
Shklovskii Anderson insulator layers (solid line). As can be 
seen, the last prediction describes the data best. Inset: T 2 - 
scaled transresistance per square for the same sample. 
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asAfarZ 1 « Q « r 1 ) ~ C 2 (e 2 /^)(^/^o)^ 2 ^ 3 (4) 

where r w = £ln(wo/w), wo = fceT/fi, and Ci, i = 1,2 
are numerical constants of order unity. Substituting 
equations ([3]) and ([4]) into ([2]) and then (QJ, we find 
that the dominant contribution to the transresistance, 
if the layers are treated as effectively 3-d for q > r^ 1 
and we can assume weak static screening, would change 
with T as: p d oc T 4 exp[(T /T) 1 / 2 ]. The observed p d oc 
T 2 exp[(To/T) 1 / 2 ] temperature dependence can however 
arise if one or both of the following are the case: 

(1) The relevant localization length is much smaller than 
that obtained from the T-dependence of the DC conduc- 
tivity so that: 

^Tq/T) 1 / 2 < d. (5) 

In this way the r^ 1 < q momentum transfers (which are 
otherwise dominant) are suppressed, and the q < r^ 1 
contribution determines the low-T transresistance. Sub- 
stituting equation (j4|) into ([2|) and fl}, the transresis- 
tance would now change with temperature as: ^ oc 
T^MiTo/T) 1 ' 2 ]. 

(2) Finite-w transport is effectively 3-d due to pair-arms 
(r u ) reaching into the barrier and the second layer. From 
the strongly localized nature of electrons in the barrier 
layer, the effective £ would be much smaller and condition 

may be satisfied, even if the effective layer separation 
also decreases significantly. 

We note that since several of the relevant length scales 
are comparable (r c ~ £ ~ W ~ 2d), a small modifica- 
tion of the effective values that these parameters take 
(due to e.g. finite-size effects, correlated hopping or sur- 
face effects) could cause a change between the £ > d and 
£ < d regime in the temperature range probed, resulting 
in a different T-dependence than that predicted. We also 
note that we do not observe the expected transition from 
the q < r^ 1 to the q > r^ 1 regime, as the temperature 
is increased. However, if the apparent change in the T- 
dependence of the transresistance (see figure |4j at T«4K 
occurs when the q > r" 1 regime kicks in, then we predict 
that £ w {AK(iK/e 2 )d 2 ~ 3A. This would give a transre- 
sistance that goes as p d oc T 2 exp[(T /T) 1 / 2 ] for T > 4K, 
or q < r~ , as we have observed. 

In conclusion, we have observed the Coulomb drag ef- 
fect between two 200 A thick insulating a-Sii-^Nb^ films 
(with 0.07 < x < 0.08) separated by a 50 - 200A thick 
SiO based barrier. We were able to retrieve accurate 
linear-response data for x = 0.070(±0.002)%, and layer 
separations of 100(±10)A. The temperature dependence 
of the transresistance in such samples was found to be 
in agreement with that predicted for ES Anderson In- 
sulator films, provided that the localization length in 



the a-Sii-zNb^ layers is smaller than that inferred from 
the temperature dependence of the DC layer-resistances. 
Our study suggests that whilst theoretically the Coulomb 
drag effect is a useful technique for distinguishing the 
insulating states of thin films, it is experimentally chal- 
lenging due to the complex dielectric properties of disor- 
dered thin films at low energies, which result in non- linear 
inter- and intra- layer excitations becoming dominant at 
practical temperatures and sample dimensions. Experi- 
mental studies of the non-linear crosstalk regime between 
thin insulating films, along with simulations of the non- 
linear (current-dependent) transresistance in such sys- 
tems, may prove to be the most productive method of 
studying the detailed nature of the observed excitations. 
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